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Abstract 

Two approaches to nonperturbative renormalization are discussed for theories 
quantized on the light cone. One is tailored specifically to a calculation of 
the dressed-electron state in quantum electrodynamics, where an invariant- 
mass cutoff is used as a regulator and a Tamm-Dancoff truncation is made 
to include no more than two photons. The other approach is based on Pauli- 
Villars regulators and is applied to Yukawa theory and a related soluble model. 
In both cases discretized light-cone quantization is used to obtain a finite 
matrix problem that can be solved nonperturbatively. 

I. INTRODUCTION 

Light-cone quantization |lj has attracted some interest as a means to perform nonper- 
turbative analyses of quantum field theories 0. There are good reasons to hope that this 
technique will provide the leverage needed to obtain a qualitative, and perhaps quantitative, 
connection between quantum chromodynamics (QCD) and the constituent quark model [Q. 
Given the complexity of QCD, it is useful to first study simpler theories such as quantum 
electrodynamics (QED) and even models in 1 + 1 spacetime dimensions rather than 3 + 1 
dimensions. 

Bound-state calculations in QCD 3+1 and QED 3+1 require nonperturbative renormaliza- 
tion. Most attempts at such calculations have used Tamm-Dancoff truncations [|J and 
cutoff-type regularization, which require counterterms that depend on Fock sector ||. An 
example of such a calculation is given here for the electron's anomalous moment ||. We then 
explore the practicality of Pauli-Villars regularization as an alternative. In particular, 
we consider a simple heavy-fermion model abstracted from the Yukawa model. 

We define light-cone coordinates by 
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x ± = t ± z , x.± = (x,y). (1.1) 
Momentum variables are similarly constructed as 

P ± = E± Pz , P± = (p X ,Py). (1.2) 

The dot product is written 

p ■ x = -(p + x~ + p~x + ) — p± • xj_ . (1.3) 

The time variable is taken to be x + , and time evolution of a system is then determined by 
the conjugate operator V . The energy E is replaced by the light-cone energy p~, and the 
momentum p by the light-cone momentum p = (p + , p±). The light-cone Hamiltonian is 

H LC = V + V- - V\ , (1.4) 

where V + and V± are momentum operators conjugate to x~ and Xj_. The eigenvalue problem 
is 

H LC ^ = M 2 ^ , V3> = , (1.5) 

where M is the mass of the state. 

Some of the advantages of light-cone coordinates are the following: They admit the 
largest possible set of nondynamical generators. In particular, boosts are kinematical. For 
many theories of massive particles, the perturbative vacuum is the physical vacuum, because 
pf = \Jp 2 + m 2 +p z > implies that no particle state can contribute to the P + = vacuum. 
Thus there is no need to compute the vacuum state before computing massive states. Also, 
well-defined Fock-state expansions exist, with no disconnected vacuum pieces. 

Such expansions are written as 

^ = E f[dx] n [d 2 k ± ] n ^ n (x } k ± )\n:xP + ,xP ± + k ± } } (1.6) 

n ^ 

with n the number of particles, % ranging between 1 and n, (P + , P±) the total light-cone 
momentum, and 

[dx] n = 4^6(1 -^ Xi )ll-^=, [d 2 k ± } n = An^Y, k±0llxr- (1-7) 

i=l i=l ^y^i i=l i=\ ^ 

In the Fock basis {\n : pf, p_u)}, V + and V± are diagonal. The amplitude ip n is interpreted 
as the wave function of the contribution from states with n particles. 

A common numerical technique is discretized light-cone quantization (DLCQ) ||, in 
which periodic boundary conditions are assigned to bosons and antiperiodic to fermions in 
a light-cone box — L < x~ < L, —L± < x,y < Lj_. Integrals are replaced by trapezoidal 
approximations on a grid: p + — > jn, pj_ — > (^-n x , j^n y ), with n even for bosons and 
odd for fermions. The limit L — > oo can be exchanged for a limit in terms of the integer 
resolution K = -P + . The longitudinal momentum fraction x = p + / P + becomes n/K. H^q 
is independent of L. 
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Because the rii are all positive, DLCQ automatically limits the number of particles to 
no more than ~K/2. The integers n x and n y range between limits associated with some 
maximum integer Nj_ fixed by L± and a cutoff that limits transverse momentum. 

To reduce the size of the discrete matrix problem, a Tamm-Dancoff truncation |3[] in 
the number of particles can be applied. This has serious implications for renormalization. 
These include severe sector dependence of counterterms ||, and, for QED, violation of the 
Ward identity. 

Regularization via cutoffs typically involves limits on the invariant mass. A limit can be 
placed on the total invariant mass of each Fock state 



m 



2 + k 2 



±i 



i 



<A 2 



(1.8) 



or on the invariant mass of each particle 

m 2 + k 2 ,. 



<A 2 . 



X j 



;i.9) 



There can also be a limit on the change in invariant mass across each matrix element of H-^q 



+ k 2 



m m2 j + k ±j 



i 



.1 .; 



E 



< A 2 



;i.io) 



II. THE ANOMALOUS MOMENT 

The anomalous moment a e = F 2 (0) can be computed from a spin-flip matrix element of 
the electromagnetic current 

-F 2 (q 2 ) = -±-(P + q,UJ + mP,l) (2.1) 



2m e va ' 2P^ 

in the standard light-cone frame q = (0, qj_/P + , qj_ = q\x). Brodsky and Drell |I(J have 
given a useful reduction of this matrix element to the form 





a e = -2m e ^]e j ^ J [dx] n [d 2 k±\ n ip* n] (x, k±) ^ TvT— i>ni(x, k±) , (2-2) 

j n i^j 1» 

where ej is the fractional charge of the struck constituent and Xi = pf/P + . The wave 
functions ip n satisfy coupled integral equations obtained from H^c^ = M 2 ^> . The QED 
light-cone Hamiltonian has been given by Tang et al. [IT]. However, the bare masses and 
couplings must be computed from sector dependent renormalization conditions. 

Consider the case where there are at most two photons and only one electron. The 
Fock-state expansion can be written schematically as 

^ = ^o|e) + ^i|e7) + ^ 2 |e77). (2.3) 
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Here ipi and ip 2 are column vectors that contain the amplitudes for individual Fock states 
with one and two photons, respectively. The eigenvalue problem becomes a coupled set of 
three integral equations 

m^o + b{ • & + • i) 2 = M 2 ^ , 
b^o + AuA + A l2 $ 2 = M% , (2.4) 
b 2 Vo + Atfli + = M 2 ^2 , 

where m is the bare electron mass and b| and Aij are integral operators obtained from 
matrix elements of H^q. 

The bare electron mass in the one-photon sector is computed from the one-loop self 
energy allowed by the two-photon states. We then require that mo be such that M 2 = ml 
is an eigenvalue. The second and third equations can be solved for ipi/ipa and "0 2 /V'o- Then 
the first equation yields uiq. Normalization of \& fixes the value of ipo- 

The bare coupling for the electron-photon three-point vertex depends on the initial and 
final momenta of the electron and on the sectors between which the coupling acts. The 
momentum dependence is present because the amount of momentum available constrains 
the extent to which higher order corrections can contribute. Similarly, the sector depen- 
dence makes itself felt when the number of additional particles in higher-order corrections 
is restricted. The coupling is fixed by the ratio of the e7 — >e transition matrix element to 
the bare vertex at zero photon momentum. 

In the present calculation we use a Tamm-Dancoff truncation to {e, ery, 677}, a nonzero 
photon mass m 7 = m e /10, and a moderate coupling a = 1/10. Some results are given 
elsewhere ||. When only states with at most one photon and no pairs are retained, one can 
show that a e reduces to 

ami r dxd 2 k ± 6(A 2 - (m e 2 + kj)/x - (mj + fcj)/(l - x)) 
it 2 J x [ml — (ml + k^/x — {ml + fc^)/(l — x)} 2 



which in the limit of A — ► 00 becomes |10| 



a r 1 2x 2 (l — x)dx 



2,ti Jo x 2 + (1 — x)(m 1 /m e ) 2 
For m 7 = 0, this yields the standard Schwinger contribution [O of a/2u. 



(2.6) 



III. YUKAWA THEORY AT ONE LOOP 

As an alternative approach to regularization, we consider Pauli-Villars [0] regularization 
of the 3 + 1 Yukawa model [0,113. The one- loop fermion self-energy is proportional to 



/(/A M 2 ) = -±. ^ t,, ini r n ; ^ - D X ) , (3.1) 



dl + d 2 l ± (g+) 2 lj + (2g+ - l + ) 2 M 2 
fi 2 J l+{q+ -l+) 2 M 2 -D 1 

where q is the fermion momentum, fi is the boson mass, M is the fermion mass, and 
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TABLE I. Values of the subtracted integral I su b(M 2 /^ 2 , fiff/j, 2 ) in the limit of infinite cutoff. 
The Pauli-Villars masses are /i 2 = 10/i 2 , /i 2 = 50/x 2 and /i 2 = 100/j, 2 . 



M 2 





0.05/i 2 


O.lfi 2 


0.2/i 2 




-0.064 


0.70 


1.37 


2.70 



_ M 2 + lj M 2 + l 2 

J^l — , i — r 



(q + -l + )/q + ' 

The boson mass fi sets the energy scale. When M 2 = we obtain 



(3.2) 



/(/i 2 ,0) 



7T 



A 2 // /A 2 ^ 



(3.3) 



In order to maintain /(/i 2 , M 2 ) oc M 2 , three Pauli-Villars bosons are needed: 

/ sub (/i 2 , m 2 , = /(/i 2 , m 2 ) + £ a/(/i 2 , M 2 ) . 



(3.4) 



i=i 



The Cj are chosen to satisfy 



(3.5) 



i=i 



i=i 



i=l 



A DLCQ calculation of 7 su b has been done [16], with values of 20, 22, and 24 for K and 
25 through 30 for A^. Modification of the trapezoidal rule, with introduction of unequal 
weights, is necessary to obtain sufficient accuracy. Each integral in ( |3.4| ) was separately 
extrapolated to infinite K and A^j_ via fits to either Co + ai/K 3 + bi/Nj_ or cq + ax/K 3 + 
a 2 /K A + bi/N]_ + b 2 /N\. The latter was used for the \i\ integral. Extrapolation after 
subtraction is not as accurate. The resulting values of J SU b were extrapolated to infinite 
cutoff by fits to a + 6/A 2 . These fully extrapolated values are given in Table |. 

The magnitude of the error in each extrapolated integral was found to be < 0.02 when 
compared to the analytic result for M 2 = 0. This implies an error of ±0.04 in the J su b values. 
The extrapolation in A 2 induces additional uncertainty reflected in the miss of zero by 0.06 
for M 2 = 0. The values in Table | are consistent with 7 su b oc M 2 to within this amount of 
error. 

The number of Fock states required for Pauli-Villars particles is approximately 1.5 times 
the number for physical states. A listing of counts for two cases is given in Table |T[ Making 
Hi larger does decrease the number of Pauli-Villars states but this increases the coefficients Cj 
and thereby amplifies errors in the integrals. Also, with fewer states, the integrals themselves 
are approximated less accurately. 

We could also consider the boson self energy. To lowest order there is a fermion loop 
contribution 



dl + d 2 l ± q + (lj + M 2 ) 6(A 2 -D 2 
ALL\ l+ 2 (q+-l+) 2 fi 2 -D 2 



(3.6) 



5 



TABLE II. Number of Fock states used in two typical cases. 



physical Pauli-Villars boson states 



A 2 In 2 


K 


N±_ 


boson states fi 2 = 10/i 2 


u 2 = 50u 2 


ul = 100u 2 


total 


200 
200 


20 
24 


25 
30 


25975 22602 
44943 39162 


11142 
19293 


3305 
5695 


37049 
64150 


where 






D 2 ^q +2 (M 2 + ll)/[l + (q + - 


/ + )], 




(3.7) 


and a < 


p 4 contribution 
















r dl + d 2 l ± dk + d 2 k ± 6 (A 2 


- D 4 ) 




(3.8) 








J q+l+k+(q+ -1+ -k+) fi 2 - 


-D A ' 





where 

n - + *1 At 2 + A;j /i 2 + (/ ± + fc ± ) 2 
4 ~ (g+-Z+-A;+)/g+' 



A Pauli-Villars fermion may be needed. 



IV. A HEAVY-FERMION MODEL 



By some severe modifications of the Yukawa Hamiltonian |T7[ we obtain the following 
model Hamiltonian: 



eff 
LC 



dq + d 2 q± 



+9 



dpid 2 p±i r dp 2 d 2 p±2 f dq + d 2 q± 



fi 2 + q 2 ± f 
— — aja. 



/jj + gj t 

g+ 1 - 



(4.1) 



x 



4<5(p 1 - p 2 + g) + a £ 5(p 1 -p 2 -q) 

+ia\ 1 5(p 1 ~P 2 + q)+ iaiqS^ -p 2 ~q) 

The kinetic energy of the fermion is no longer momentum dependent and only a modified no- 
flip three-point vertex remains as an interaction. The fermion then acts as a "static" source 
for the boson. We include one Pauli-Villars field, which will prove sufficient in this case. 
Similar Hamiltonians, without the Pauli-Villars field, have been considered in equal-time 
|T%| and light-cone coordinates [fH)| . 

We write the eigenvector as a Fock-state expansion 



n,ni 



^16n s p 



8=1 



d q?d 2 q ±i r dr+d 2 r ±J 



(4.2) 



' (n,ni) 



(gj ' ^ 5 E) -y^f 6 £<x II 4, IT 4 . 1 0) 
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normalized according to $J • ^ a = 16n 3 P + 5(P' — P), which yields 



I (n,ni) 



(^,?: i ;p-£g.-£ 5 



n „ m . 

1 = e n / n / ^+rfv ±J 

n,m i J j J 

For $0- to satisfy the Schrodinger equation (|1.5| ), the amplitudes must satisfy 

dn,ni) 



(4.3) 



M 2 -M 2 _£^±^i_£^ + r ^ 



y% 



(4.4) 



+ ^E 



(n— l,ni) 



'I67r 3 g+ 



Vl67r J r 
1 



n 



167r 3 r^ 



(n,ni — 1) 



) — J — 1 5 — J + l ) ' ' ' ) 



Ln V P) 



The structure of this coupled set of integral equations is deliberately identical in basic form 
to the equations considered by Greenberg and Schweber |TJ|] . Therefore, we transcribe their 
ansatz for a solution to light-cone form 



I (n,ni) 



r ff) n Hfl) ni -Q 



9i 



n 



This does work as a solution if Mq(//i) is chosen to satisfy 



(4.5) 



M 2 - M 2 



167T 3 



dyd 2 q± f dzd 2 r± 



From the normalization condition (|4.3|) we obtain 



i = exp 



167T" 



ydyd 2 q_ L 



+ 



zdzd 2 r± 
(/if + r\f 



(4.6) 



(4.7) 



The bare mass and wave function renormalization are thus determined as functions of the 
Pauli-Villars mass. 

To fix the coupling we could use the slope of the fermion no-flip form factor, which is 
related to the transverse size of the dressed fermion. The form factor is most easily evaluated 
from [ID | 



F(Q 2 



1 



2P+ 



(p + p 7 T |J + (o)|Pt) 



(4.8) 



EeJ 16^5(1 -£*<)*(£ Mil 



dxid 2 p±i 
16tt 3 



where the matrix element has been evaluated in the frame with 



M 



P = (P + , P~ = — , ± ) , p, = (0,p~ = 2p, ■ P/P + , p 7± ) , Q 2 = p 2 ± , 
ej is the charge of the jth constituent, and 



P±i - Xip 1± i ^ j 
P±i + (1 - Xi)p 1± i = j . 



(4.9) 



(4.10) 



A sum over Fock states is understood. 

When the fermion is assigned a charge of 1, and the bosons remain neutral, the analytic 
solution for the amplitudes yields 

^-H'/^&A^M- (4 ' u) 



with 



From this we find 



q± = q± - yp~f± ■ 



F'(0) = -9 2 j 



dyd 2 q± y 3 



2/i 2 



fi 2 + ql 



+ P-V term . 



16tt 3 (fi 2 + q 2 ± ) 3 

Numerically, the slope is computed from a finite-difference approximation to 



F'(0) = E II / dqtd 2 qu\{ / dr+d 2 r ±J 



x 



z-_ 2 
4 - ±j 



(n,ni) 



(^■;£-E^-E^ 



X 



(n,m). 



%,rj,p- E^-E^) 



(4.12) 



(4.13) 



(4.14) 



With the bare parameters determined, we "predict" a value for (:0 2 (O):). For the analytic 
solution, this expectation value reduces to 



<:0 2 (O):) 



if 



8ttV 



i-^-^V 2 



(4.15) 



A 2 A 2 A 2 _ 

From a numerical solution it can be computed from a sum similar to the normalization sum 

Ell / dqtd 2 qu II / ^V ± , f E - 

n=l,ni=0 j jr 



(:0 2 (O):) 



x 



{n ' ni) (^i,;Z-E^-E^ 



(4.16) 
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V. SUMMARY 



For the anomalous moment calculation there remain several hurdles. The Tamm-Dancoff 
truncation results in logarithmically divergent four-point graphs. To deal with these will 



probably require use of scattering processes, such as Compton scattering p0| , to obtain 
renormalization conditions. Verification of the removal of all logarithms and restoration of 
symmetries can then be undertaken. Also neglected up to this point have been zero modes, 
photon modes of zero longitudinal momentum |[21||. How they might be included has been 



indicated by Kalloniatis and Robertson [22 



Additional physics could be included in the calculation by introducing an effective in- 
teraction from Z graphs or even putting eee + states in the basis. In the latter case, photon 
mass renormalization must be done. 

In the Yukawa-model calculations we have learned that Pauli-Villars Fock states increase 
the basis size by only 150%, which may not be prohibitive. To perform such calculations ac- 
curately with a minimal basis size, improvement of ordinary DLCQ, by inclusion of weighting 
factors, is critical. 

We have found a simple 3 + 1 model, related to Yukawa theory, which can be solved 
analytically. Here we will attempt a nonperturbative numerical solution to further test the 
use of Pauli-Villars regularization in DLCQ. If successful, we can begin to increase the 
complexity of the model, eventually reaching the full Yukawa theory. 
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